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■■ Abstract 



A three-dimensional, space-frequency model for the excitation of electromagnetic 



> 

^ . radiation in a free-electron laser is presented. The approach is applied in a numerical 

^ ■ particle code WB3D, simulating the interaction of a free-electron laser operating in 

TJ" I the linear and non-linear regimes. Solution of the electromagnetic excitation equa- 



tions in the frequency domain inherently takes into account dispersive effects arising 
from the cavity and the gain medium. Moreover, it facilitates the consideration of 
^i statistical features of the electron beam and the excited radiation, necessary for the 

Oh' study of spontaneous emission, synchrotron amplified spontaneous emission (SASE), 

super-radiance and noise. 

;h ' We employ the code to study the statistical and spectral characteristics of the 

radiation generated in a pulsed beam free-electron laser operating in the millimeter 
wavelengths. The evolution of radiation spectrum, excited when a Gaussian shaped 
bunch with a random distribution of electrons is passing through the wiggler, was 
investigated. Numerical results of spontaneous emission power along the wiggler are 
compared to analytical predictions in the linear regime. In the first few periods, 
the power is excited from shot-noise in the low-gain regime. An exponential growth 
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of SASE in the high-gain regime is inspected after passing a sufficient number of 
periods, until saturation occurs when arriving to the non-Hnear regime of the FEL 
operation. 



1 Introduction 

Electron devices such as microwave tubes and free-electron lasers (FELs) uti- 
lize distributed interaction between an electron beam and electromagnetic ra- 
diation. Random electron distribution in the e-beam causes fluctuations in cur- 
rent density, identified as shot noise in the beam current [l]-[4]. Electrons pass- 
ing through a magnetic undulator emit a partially coherent radiation which is 
called undulator synchrotron radiation [5]. The electromagnetic fields excited 
by each electron add incoherently, resulting in spontaneous emission [6]-[13]. 
When the electron beam is modulated or pre-bunched, the fields excited by 
electrons become correlated, and coherent summation of radiation fields from 
individual particles occurs. If all electrons radiate in phase with each other, 
the generated radiation becomes coherent {super-radiant emission). 

In high-gain FELs, utilizing sufficiently long undulators, the spontaneous emis- 
sion radiation excited in the first part of the undulator is amplified along 
the reminder of the interaction region resulting in self-amplified spontaneous 
emission (SASE) [14]-[18]. Super-radiant emission emerges if the electrons are 
injected into the undulator in a single short bunch (shorter than the oscilla- 
tion period of the emitted radiation) [19]-[25] or enter as a periodic train of 
bunches at the frequency of the emitted radiation [26]-[30]. 

Investigation of spontaneous and super-radiant emissions, as well as SASE, 
call for analytical and numerical models, that can describe non-stationary 
stochastic processes involved in the generation of incoherent or partially co- 
herent radiation. In addition to the statistical features of the electron beam 
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and radiation, the model should take into account dispersive effects evolving 
in the medium, in which the radiation is excited. The model presented in this 
paper, utilizes an expansion of the total electromagnetic field (radiation and 
space-charge waves) in terms of transverse eigenmodes of the medium, in which 
the field is excited and propagates [31,32]. The interaction between the elec- 
tron beam and the electromagnetic field is described by a set of coupled-mode 
excitation equations, that expresses the evolution of mode amplitudes in the 
frequency domain. The excitation equations are solved self-consistently with 
the equation of particles motion, which describes the electron beam dynamics. 



2 Presentation of the electromagnetic field in the frequency do- 
main 

The electromagnetic field in the time domain is described by the space-time 

electric E(r, t) and magnetic H(r,t) signal vectors, r stands for the {x,y,z) 

coordinates, where {x, y) are the transverse coordinates and z is the axis of 

propagation. The Fourier transform of the electric field is defined by: 

+00 
E(r,/)= J E{r, t) e+^'-f' df (1) 

—00 

where / denotes the frequency. Similar expression is defined for the Fourier 
transform H(r, /) of the magnetic field. Since the electromagnetic signal is real 

(i.e. E*(r,t) = E(r,t)), its Fourier transform satisfies E*(r,/) = E(r, — /). 

Analytic representation of the signal is given by the complex expression: 

E(r,t) = E(r,t)-jE(7;t) (2) 



where 



+00 , ., 
E^t) = / ^^ dt' (3) 



-00 



is the Hilbert transform of E(r, t). Fourier transformation of the analytic rep- 
resentaion (2) results in a 'phasor-like' function E(r, /) defined in the positive 
frequency domain and related to the Fourier transform by: 



E(r,/) = 2E(r,/).ii(/) 



2E(r,/) />0 

(4) 
/<0 



The Fourier transform can be decomposed in terms of the 'phasor like' func- 
tions according to: 

E(r,/) = iE(r,/) + iE*(r,-/) (5) 

and the inverse Fourier transform is then: 

+00 f oo \ 

E(r, t)= J E(r, /) e-^'^^' df = ^U E(r, /) e-^'^^' df (6) 

-oo lo J 



3 The Wiener-Khinchine and Parseval theorems for electromag- 
netic fields 

The cross-correlation function of the time dependent electric E(r,t) and mag- 
netic H(r,t) fields is given by: 

+ 00 

Reh{z, r) = f I f f [E(r, t + r) x H(r, t)] ■ z dx dy\ dt (7) 



-oo 



Note that for finite energy signals, the total energy carried by the electromag- 
netic field is given by W{z) = Reh{z, 0) 

According to the Wiener-Khinchine theorem, the spectral density function of 
the electromagnetic signal energy Seh{z, f) is related to the Fourier transform 
of the cross-correlation function Reh{z, t) through the Fourier transformation: 
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+ 00 

Seh{z, f)= J Reh{z, r) e+^^^f^ '^^ = J J t^^^' /) ^ H*(^' /)] -^dxdy 



JJl[E{r,f)xH*{r,f)]-zdxdy />0 

// i [E(r, -/) X H*(r, -/)!* ■ z dxrfy / < 



(8) 



Following Parseval theorem, the total energy carried by the electromagnetic 
field can also be calculated by integrating the spectral density Seh{z, f) over 
the entire frequency domain: 



W{z) = J 

— oo 


Seh{z, f) df 


oo 

-I 




[II>{^^ 


identify: 


d\ 


""^^^-'-^l 



df 



(9) 



df 2 



E(r,/)xH*(r,/) -zdxdy 



(10) 



as the spectral energy distribution of the electromagnetic field (over positive 
frequencies) . 



4 Modal presentation of electromagnetic field in the frequency do- 
main 

The 'phasor like' quantities defined in (4) can be expanded in terms of trans- 
verse eigenmodes of the medium in which the field is excited and propagates 
[32]. The perpendicular component of the electric and magnetic fields are given 
in any cross-section as a linear superposition of a complete set of transverse 
eigenmodes: 



Ex(r,f) = 5:[C+,(z,/)e+^-'=-^ + C^,{z,f)e-^'^^^] S,^{x,y) 



H^(r,f) = ^[C+,(z,/)e+^-^-^ - C^,{z,f)e-^'^^^] n,^{x,y) (11) 



C+q{z, /) and C_q{z, f) are scalar amplitudes of the qth forward and backward 
modes respectively with electric field Sqj_{x,y) and magnetic field 'Hq±{x,y) 
profiles and axial wavenumber: 



Kzq — S 



jjkj_q — k"^ k < k±q (cut-off modes) 



(12) 



k'^ ~ k±q k > k±q (propagating modes) 



Expressions for the longitudinal component of the electric and magnetic fields 
are obtained after substituting the modal representation (11) of the fields 
into Maxwell's equations, where source of electric current density J(r, /) is 
introduced: 



E,(r, /) = ^ [C^q{z, f) e+^-'^-^ - C^q{z, f) e-^-'^-^l ^(x, y) + — -- J.(r, /) 
H.{r, f) = J2 [C+qiz, f) e+^-'=-^ + C_,(z, /) g-^'^-^] 7?,,(x, y) (13) 



The evolution of the amplitudes of the excited modes is described by a set of 
coupled differential equations: 



dz 



C±q{zJ) 



■T 



2K 



pTj'^zq^ 



^ J^(r,/)±zJ,(r,/) 



S*{x,y)dxdy (14) 



The normalization of the field amplitudes of each mode is made via each 
mode's complex Poynting vector power: 



K 



£q±{x,y) y<'H*q^{x,y) -zdxdy 



(15) 



and the mode impedance is given by: 



V ^0 


k 


/mo 

V ^0 


k 



fiza 



27r/eo 



for TE modes 
for TM modes 



(16) 



Substituting the expansion (11) in (10) results in an expression for the spectral 
energy distribution of the electromagnetic field (over positive frequencies) as 
a sum of energy spectrum of the excited modes: 



dWjz) 
df 



I 



E o \C\,izJ)\' - \C^,izJ)\' ^{Af,} + 



Propagating 



+ E Q{C^,iz,f)Cl^iz,f)}Q{Ar,} (17) 

Cnt-o// 

The power spectral density carried by the propagating mode q during a tem- 
poral period T is given by the ensemble average: 



^ = ^i^.(-,/)rn 



(18) 



where Pg = iK{Arj. 



5 The electron beam dynamics 



The state of the particle i is described by a six-components vector, which 
consists of the particle's position coordinates Tj = {xi,yi,Zi) and velocity vec- 
tor Vj. The velocity of each particle, in the presence of electric E(r,t) and 
magnetic B(r,t) = /ioH(r,t) fields, is found from the Lorentz force equation: 



dvi 1 f e 1 ,^, , ^, ,. d'ji 



(19) 



where e and m are the electron charge and mass respectively. The fields repre- 
sent the total (DC and AC) forces operating on the particle, and include also 
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the self-field due to space-charge. The Lorentz relativistic factor 7^ of each 
particle is found from the equation for kinetic energy: 

d'li e 1 ^, , , , 

-7^ = ^— v,-Er„t 20 

where c is the velocity of light. 

The time it takes a particle to arrive at a position z, is a function of the 
time toj when the particle entered at 2; = 0, and its instantaneous longitudinal 
velocity Vzi{z) along the path of motion: 

'^ ' 

The current distribution is determined by the position and the velocity of the 
particles in the beam: 

J(r, t) = -g ^ [^ (5(x - X,) 5{y - y^) 5 [t - U{z)] (22) 

here q = ^ is the charge of each of the N macro particles in the simulation 
(Jo is the DC current of the e-beam pulse of temporal duration T). The 'phasor 
like' current density is given by: 



+00 
3{r,f) = 2u{f) J 3{r,t)e-^'''f'dt = 

= -2g E (—] ^(^ - ^0 S{y - Vi) e-^-2-^*»(^) u{f) (23) 



1=1 \ ^I , 



6 Numerical results 

The WB3D code was used to investigate the excitation of spontaneous emis- 
sion in a millimeter wave free-electron maser (FEM), with operational pa- 
rameters given in Table 1. The corresponding dispersion curves of the FEM 



are shown in Fig. 1. When the beam energy is set to 1.375 MeV, there are 
two separated intersection points between the beam and waveguide dispersion 
curves, corresponding to the "slow" {Vg^ < Vzg) and "fast" {Vg^ > Vz^) synchro- 
nism frequencies 29 GHz and 100 GHz, respectively. Lowering the beam en- 
ergy to 1.066 MeV, results in a single intersection at 44 GHz ("grazing limit"), 
where the beam dispersion line is tangential to the waveguide dispersion curve 

The evolution of spontaneous emission power spectrum in the vicinity of the 
upper synchronism frequency 100 GHz is drawn in Fig. 2. a. The power growth 
along the wiggler as a function of the wiggling periods A^^, is described in Fig. 
2.b. In the first few periods, the mutual interaction between the electromag- 
netic radiation and the electron beam is small and the power amplification 
is low. Within this stage, the spontaneous radiation power increases propor- 
tional to A^^. An exponential growth of SASE is inspected later after passing 
a sufficient number of periods, revealing that the interaction enters to the high 
gain regime, until saturation occurs when arriving to the non-linear regime of 
the FEL operation. Fig. 3 describes the power evolution in the case of grazing. 
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Table 1 

Operational parameters of millimeter wave free-electron maser. 

Accelerator 

Electron beam energy: £"^ = 1^3 MeV 

Electron beam current: /o = l A 

Wiggler 

Magnetic induction: 5^=2000 G 

Period: A^=4.444 cm 

Waveguide 

Rectangular waveguide: 1.01 cm x 0.9005 cm 

Mode: TEqi 
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Figure 1. FEM dispersion curves. 
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Figure 2. Evolution of spontaneous emission: (a) spectrum and (b) power along the 
wiggler. The beam energy is £"^=1.375 MeV. 
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